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We study the sensitivity and resolution of phase measurement in a Mach-Zehnder interferometer 
with two-mode squeezed vacuum (n photons on average). We show that super-resolution and sub- 
Heisenberg sensitivity is obtained with parity detection. In particular, in our setup, dependence of 
the signal on the phase evolves n times faster than in traditional schemes, and uncertainty in the 
phase estimation is better than 1/n. 
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Different physical mechanisms contribute to phase 
measurement. Thus, measuring phase provides insight 
into a number of physical processes. Therefore, improved 
phase estimation benefits multiple areas of scientific re- 
search, such as quantum metrology, imaging, sensing, 
and information processing. Consequently, enormous ef- 
forts have been devoted to improve the resolution and 
sensitivity of interferometers. Sensitivity is a measure of 
the uncertainty in the phase estimation, while resolution 
is rate at which signal changes with changing phase. 

In what follows, we direct our attention to quantum 
interferometry. The benchmark that quantum interfer- 
ometry is compared against is one with coherent light 
input and intensity difference measurement at the out- 
put of a Mach-Zehnder interferometer (MZI) . In general, 
phase sensitivity of this benchmark is shot-noise limited, 
namely A(p = where n is the average number of 

photons. However, better sensitivity is possible if nonlin- 
ear interaction between photons in the MZI takes place 
[l[ • In what follows, we only consider phase accumulation 
due to linear processes. 

In 1981, Caves pointed out that by using coherent light 
and squeezed vacuum one could beat the shot-noise limit 
Aip < (super-sensitivity) Q. In the work of Boto 

et ai, it was shown that by exploiting quantum states 
of light, such as NOON states, it is possible to beat the 
Rayleigh diffraction limit in imaging and lithography (su- 
per resolution) while also beating the shot-noise limit in 
phase estimation 

(ESSS- Finally, it was shown in Ref. 
[7| that input state entanglement is important in order 
to achieve super-sensitivity in a linear interferometer. 

non-classical light Experimental realization of these 
predictions have been hindered by the fact that en- 
tangled states of light, with large numbers of photons, 
are difficult to obtain. Therefore we turn our atten- 
tion to the brightest (experimentally available) non- 
classical light — two- mode squeezed vacuum (TMSV). 
A state of TMSV is a superposition of twin Fock states 
l^n) = J2^=o VPn (") \ n , n ), where the probability of 
a twin Fock state \n,n) = \n) A \n) B to be present de- 




FIG. 1: Mach-Zehnder interferometer used in the calcula- 
tions. Two-mode squeezed vacuum input state \ipn) is in- 
dicated together with the intermediate state \ip} = PipU \ ipn)- 
Vertical dash-dot lines indicate places where two measure- 
ments fiAB and YIa are to be implemented. 



pends on average number of photons in both modes of 
TMSV, n, in the following way p n (n) = (1 — t n )t^ with 
tft = l/(l + 2/n) !■ 

Light entering a MZI in TMSV state exits a loss- 
less interferometer in the state \ipf) = Umzi \ipn), where 
the MZI is described by the unitary transformation 
Umzi (Fig. [1]). This transformation, in terms of 
the field operators for the optical modes a and b, 
is Umzi = UP V U = exp (tp (a^b — S^aJ /2^, where 

P v —e-xjp(^—i(pGj describes accumulation of a phase dif- 
ference <p; and U =exp(ij (a^b + describes the 50- 

50 beam splitter, with a ir/2 phase shift for the reflected 
mode. In a linear medium the generator of phase evolu- 
tion is G — [fiA — Hb) /2, where ua — o^a- and fiB = b^b 
are the photon number operators in each mode. 

Phase estimation is based on the detection of light at 
the outputs of MZI. Not all detection schemes are capable 
of exploiting the full potential of non-classical light to be 
super-sensitive and super-resolving. For example, inten- 
sity difference measurement, which is standard for optical 
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interferometry with coherent light, is not phase sensitive 
at all if TMSV input is used [§] . In our work, we consider 
parity detection for our measuring scheme. The parity 
operator on output mode A is 11^ = exp(z7m^). Parity 
was originally discussed in the context of trapped ions by 
Bollinger et at and later adopted for optical interfer- 
ometry by Gerry and Campos [ll|, [II] • Super-sensitivity 
with this detection strategy has been shown for several 
classes of input states [3] • Finally, parity detection was 
also shown to allow better than classical resolution with 
coherent light at the shot-noise limit (SNL) 

Parity measurement on mode A at the output of MZI 
is computed by (jIa^ — (ipt\ ^Ia \"*P()- It has been shown 
in Ref. [bj ] that the parity measurement on mode A after 
the final beam splitter is equivalent to the measurement 
that is constructed from all the \M, M') -> \M', M) pro- 



jectors as follows [14 1 
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acting on the inner modes of MZI, such that \ip) = 
PipU\tpn). Our use of the (iab operator here highlights 
the fact that parity detection combined with 50-50 beam 
splitter provides a measurement scheme that includes all 
of the phase-carrying off-diagonal terms in the two-mode 
density matrix 14]. Calculation of {(lab) simplifies sig- 
nificantly once it is noted that such an operator, as well 
as a beam splitter, does not change the total number of 
photons in the state. Thus a lossless MZI with the par- 
ity detection scheme does not mix twin Fock states with 
different number of photons giving: 



n=0 



(2) 



where {J^-aJ = (n, rc|C/j[flzi-ilA^MZi \n,n) is the expecta- 
tion value of the parity operator for twin Fock state input . 
In turn, the expression (j^-A^ = ( — l) n Pn (cos (2cp)), 
given in terms of Legendre polynomials P n , could be 
found in Ref. fl~5j ] . Finally, one can identify our expres- 
sion in Eq. @ with the generating function for Legendre 
polynomials [l6| and arrive to the following: 
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ABU = 



1 + n(n + 2) sin 2 <p 



(3) 



where an additional tt/2 phase shift was introduced. 
Eq. (J3j) is the central result of this paper and, in what 
follows, it will be used to study the resolution and sensi- 
tivity of our proposed scheme. 

Let us compare here the signal outcomes of the TMSV 
scheme with n = 10 to coherent-state-based optical in- 
terferometry with n — 100 (see Fig. [2]). Intensity 




FIG. 2: Measured signals at the output of MZI with coherent 
light (n — 100) and TMSV (n = 10) inputs against accumu- 
lated phase difference. Solid and dashed lines are the outputs 
of parity measurement for TMSV and coherent light, respec- 
tively. Dotted line, given for comparison, is a scaled down 
output of intensity difference measurement on the output of 
MZI fed with coherent light. TMSV with much smaller pho- 
ton number performs as well as coherent light. 



difference measurement, with coherent state at the in- 
put of MZI, exhibits classical interference — a sinu- 
soidal dependence on the phase with an intensity inde- 
pendent period of 2-7T. In the case of parity detection 
with coherent state input, it was shown in Ref. [14] that 

11,4 ) =exp(— 2n sin 2 (v/2)) with a 2tt period and a fea- 



=exp( 

ture at the phase origin that is narrower than the classical 
curve by a factor of Sip = \ j\ffi. In our case, the width 
of the feature is further reduced by y/n + 2 times. There- 
fore, the peak in Fig. [2] is as narrow for a n = 10 TMSV 
as for a n = 100 coherent state input. 

The other aim of optical interferometry is to mini- 
mize uncertainty, A<p, of the measured phase. The low- 
est bound on the uncertainty is inversely proportional to 
quantum Fisher information Fq fl7| . In the case of a 
pure state Fq = 4 AG 2 , where AG 2 is a variance of the 
phase evolution generator given above. Thus, 
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4AG 2 ' 



(4) 



which depends on the state of the light used but not on 
the measurement. It is in this spirit that we replaced 
coherent light with two-mode squeezed vacuum in order 
to beat the sensitivity obtained in Ref. [3]. Our anal- 
ysis shows that Eq. Jl]), in the case of coherent light, 
limits the attainable sensitivity to Aipf nin — shot 
noise, while for TMSV it sets much lower limit A^ in = 
(n(n + 2)y 1 < fir 2 . This means that TMSV has a po- 
tential for super sensitive phase estimation, which has 
phase uncertainty better than 1/n and that is thus sub- 
Heisenberg. However, it remains to be seen whether sub- 
Heisenberg sensitivity could be obtained with a particu- 
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FIG. 3: Sensitivity of phase estimation obtained with parity 
measurement at <p = (dashed) against average total photon 
number. Dotes are sensitivity estimation based on quantum 
Fisher information for integer values of n. Shaded area is 
between dotted lines 1/n and l/V^. Solid line is for the 
Hofmann limit discussed in the text. Inset: Sensitivity against 
actual values of accumulated phase difference. Solid lines for 
TMSV with n = 5 and n = 25; dashed line for coherent light 
with ri = 25. 



lar measurement, namely parity. 

The variance of the phase estimation based on the out- 
come of the parity measurement can be estimated as 

= ,!,: {flAB !\ 2 , (5) 



{d (fL AB ) /d(p) 

which is a ratio of detection noise to the rate at which 
signal changes as function of phase. We have shown that 
the rate of the signal change is much higher than in the 
case of coherent state input. Thus, if parity measurement 
on the squeezed vacuum is no noisier than on the coherent 
state, sensitivity improvement is expected. 

The sensitivity of the phase estimation in our scheme 
can be estimated based on Eq. ([5]) combined with Eq. ([3]): 



A<p = 



1 + n(n + 2) sin 2 <p 



(6) 



| cos ip\ \Jn{n + 2) 

which is presented in Fig. [3] for the case of ip = 0. It 
is clear that in this case sub-Heisenberg sensitivity is ob- 
tained and that the lower bound defined by the quan- 
tum Fisher information is actually reached. Expression 
in Eq. © gives dependence of the phase sensitivity of 
our scheme on the actual phase difference as well: 
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1 + (2n(n + 2) + 1) ^ 



(7) 



y/n{n + 2) 

where expansion near the phase origin was made. This 
dependence is, in turn, compared to the one for a co- 
herent state input, which has the following functional 
dependence on the phase in the vicinity of phase origin 
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1 + (2n + 1) 



(8) 



Dependence of the phase sensitivity in both cases is pre- 
sented in the inset of Fig. [3] for n = 5 and n = 25 TMSV 
and for n = 25 coherent state inputs. Comparison shows 
that our scheme has superior sensitivity in the vicinity of 
phase origin but degrades rapidly as actual phase differ- 
ence deviates from zero. 

Better than 1/n phase sensitivity, demonstrated here 
by a linear MZI with parity detection and TMSV, is in 
violation of the so-called Heisenberg limit (HL). In what 
follows, we will argue why HL is not the true limit, if HL 
is defined as Alphl = 1/n. 

The traditional argument for the limit on the sensitiv- 
ity of the phase measurement comes from the Heisenberg 
uncertainty principle for the phase and photon number 
AipAn > 1. This argument is usually combined with 
intuitive thinking that in the case of finite energy An 
should be bounded by n, that is n > An. This argument 
is valid as long as a fixed photon number is assumed 
n = N, but has to be treated with caution for fluctuat- 
ing photon numbers. Such notions about the Heisenberg 
limit can be traced back to, for example, work by Ou 
[l|| where he speculates that the fundamental limit set 
by quantum mechanics on sensitivity is the Heisenberg 
limit Aip > 1/n since all analysis up until now had not 
shown better than 1/n sensitivity. Ou's conjecture is in 
fact not supported by his cited 1986 paper by Yurke et 
al, on SU(1,1) interferometers [3], where Eq. (9.31) was 
published with a typographical error — there is a minus 
sign where there should be a plus. 

Recent progress in quantum metrology has demon- 
strated the importance of photon number fluctuations 
for phase estimation. In order to better account for 
photon number fluctuations, Hofmann in Ref. [2(| sug- 
gested a more direct definition of the ultimate quantum 
limit of phase sensitivity A<p 2 > l/(n 2 ), where (n 2 ) 
indicates averaging over the squared photon numbers. 
Thus, in the case of high photon number fluctuations, 



An 2 



> 0, the Hofmann limit allows for 



better sensitivity of the phase measurement than the 
Heisenberg limit. Clearly (n 2 ) contains direct informa- 
tion about fluctuations where (n) 2 does not. 

In the case of our TMSV with parity, sensitivity of 
phase estimation is better than allowed by the Heisenberg 
limit, although marginally, but it is never better than 
1 /-\/2n 2 + 2n, which is the Hofmann limit. It is also never 
below the quantum Cramer-Rao lower bound set by the 
quantum Fisher information of the state at an area of 
phase accumulation. 

In order to demonstrate that the maximal phase 
sensitivity could be underestimated by the Heisenberg 
limit if photon number fluctuations are neglected, con- 
sider the following state p(n,6) = sin 2 6* |0, 0) (0, 0| + 
cos 2 6> 1 7i, n) (n, n|, which has ft = 2ncos 2 #. This state 
could appear in the context of a probabilistic twin Fock 
state generation with parity detection on a single output 
mode, since such a detection would not distinguish vac- 
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demonstration that optical interferometry with two- 
mode squeezed vacuum and parity detection provides 
an interferometric metrology strategy with sensitivity 
A(p < fiT x and resolution n^ 1 times better than the res- 
olution of standard (classical) interference. 
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FIG. 4: Phase estimation sensitivity A<p for the state p(n, 6) 
with n = 2 (top four) and n — 5 (bottom four) against 9. Solid 
lines represent sensitivity with parity measurement which sat- 
urates limit set by quantum Fisher information (dots for a 
few selected values of 8). Dashed lines represent HL sensitiv- 
ity estimation based on averaged total photon number: 1/n. 
Dot-dashed lines give Hofmann's estimation of sensitivity dis- 
cussed in the text. 



uum contribution from the twin Fock contribution when 
all photons were routed out in the other port. 

Based on quantum Fisher information, this 
state is capable of providing sensitivity A(p 2 = 
1/ (2n (n + 1) cos 2 6*), which is obtainable by parity 
measurement. This dependence of the phase sensitivity 
is presented in Fig. [4] for n = 2 and n = 5, where the 
presence of the vacuum, 9 > 0, degrades the sensitivity 
but allows for 1/n Heisenberg limit to be beat! However, 
the Hofmann limit A<p 2 = 1/ (4n 2 cos 2 0) tracks the 
phase sensitivity well; without being beaten! 

There does exist another limit based on the highest 
number of photons in the state — 1/iV, with N — 2n for 
the state considered here. However, it is not as useful as 
the Hofmann limit for a number of reasons: a) it overesti- 
mates the sensitivity as it does for p(n, 9); b) information 
about N is not readily available in experiments; c) for 
states, such as coherent and squeezed vacuum, N = oo. 

Finally, implementation of parity detection needs to be 
discussed. In proof of principle experiments, a highly ef- 
ficient photon number-resolving detector could be used. 
Such detectors with 95% efficiency and number resolv- 
ing capabilities in the tens of photons have been demon- 
strated 
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J). However, for more practical 
applications, knowledge about exact photon numbers is 
excessive. We conjecture that a scheme, which does 
not require photon counting, exists, perhaps through the 
exploitation of optical non-linearities [25j |. or projective 
measurements, and this is an area of ongoing research. 
In conclusion, the main result of this paper is our 
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